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A generalization of the max-plus transformation, which is known as a method to derive 
cellular automata from integrable equations, is proposed for complex numbers. Operation 
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rules for this transformation is also studied for general number of complex variables. As an 
application, the max-plus transformation is applied to the discrete Fourier transformation. 
Stretched coordinates are introduced to obtain the max-plus transformation whose imaginary 
part coinsides with a phase of the discrete Fourier transformation. 
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^0 , 1. Introduction 

in ■ 

£^ ■ The importance of studies on discrete systems generated from integrable equations has 



recently attracted considerable attention. Because these discrete systems conserve good prop- 
erties of integrable equations, such systems, known as integrable difference equations and 
soliton callular automata, 1 are important as models to describe engineering problems. 

Ultradiscretization is introduced as one of the systematic methods to derive discrete sys- 
tems from differential equations. 2 ' 3 It is based on an exponential transformation of field vari- 
^ \ ables u into U, 

C3\ u = e u/e , e>0, (1.1a) 

followed by the limiting process e — > +0. By virtue of a relation 

lim elog(e A/£ + e B/£ ) = max(A,B), (1.1b) 

where 

(a (A>B), 

max(A£) = < (Flc) 
[B (A<B), 

we can transfer summation and product of positive quantities a and b. Let us set a = e A / £ 
and b = e B l e respectively, and then consider the logarithms of given equations, multiply e 
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and take limit of e — > +0. Variables A and B may depend on parameter e, and if they remain 
finite at e — > +0, we have 

a + b = e A/£ + e B/£ ^max(A,B), (1.2a) 

ab = ^ A+B )/ e ^ A + B. (1.2b) 

This procedure is usually called "max-plus transformation". Many nonlinear equations have 
been transferred to full-discrete systems by this method, and the relations between them have 
been studied intensively. ' 5 

However, it is not sufficient to consider the max-plus algebra only for real variables, because 
all the variables sould be positive in order to ensure the relevance of the transformation and 
limit defined in eqs. (|1.1|) . Since the physical variables can be negative, this restriction is fatal 
from the viewpoint to discretize solutions of difference equations directly. 

In this paper, we aim to introduce a well-defined limiting process for complex variables, 
which corresponds to a generalization of the max-plus transformation for real numbers. By 
considering eqs. (|l,lj) for complex variables, we can obtain a new max-plus type algebra. We 
can discretize the transformation U of the solutions of difference equations automatically, if the 
initial values are properly chosen discrete. This method eliminates the difficulty of the max- 
plus transformation for real variables. Moreover, we can apply the max-plus transformation 
to discrete Fourier transformation (hereafter we shall abbreviate as DFT), by which we can 
solve various differential equations. 

This paper is organized as follows. In the next section, we shall introduce a max-plus 
algebra for two complex variables. A generalization of this algebla for more than two variables 
is given in §3. In §4, we apply this algebra to DFT. The final section is devoted to summary 
and discussions. 

2. Max-Plus Type Transformation for Two Complex Variables 

Let us consider the generalization of max-plus transformation for the sum of two complex 
variables u\ and u 2 '- 

u = ui + u 2 - (2.1) 

Parallel to the method for real variables, we introduce new complex quantities U\ and U 2 
corresponding respectively 

Ul = e u ^l £ and u 2 = e u ^ s , (2.2) 

where e is a positive number. Hereafter, we shall write elogu as W(e) for simplicity, and are 
going to derive its explicit expression. Substituting the relations (|2.2|) into ()2.1|) . we have 

W(e) = e log (e Ul/e + z U2,£ ) . (2.3) 
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By the definition of the logarithms on the complex plane, we can write down the real and the 
imaginary parts of the right-hand side of eq. 1)2.3)1 to see 



W(e) =eln 



e U 1 /e +e U 2 /e 



+ iearg[e Ul/e + e U2/£ 



(2.4) 



We express U\ and U 2 as Uj = Xj+iyj (j = 1, 2), where Xj, yj are real numbers, and we assume 
x\ > X2 for the nonce. First, we consider the real part of W(e). From direct calculation, we 
have 

1/2 



JJi/e + e U 2 /e 



e 2 X1 /a + e 2x 2 /e\ (1 + A) 



xi - x 2 yi-2/2 
A = seen cos . 



(2.5) 



e e 

From eqs. 1)2.4)) and 1)2.5)1 . we find the real part of W(e) as 



ReW(e) = - In 



} 2asi/e _|_ e 2x 2 /e 



(1 + A) 



= £ In (e 2xi / £ + e 2:C2 / e ) + | ln(l + A) . (2.6) 

Due to the fact x\ > x 2 and the similar calculation to the max-plus transformation of real 
variables, the first term of eq. 1)2.6)1 yields 



| In (e 2xi/£ + e 2x2/£ ) -► s x> as e -» +0. 



We note that A tends to zero as e — > 0, as long as x\ ^ x 2 is satisfied, and ln(l + A) is kept 
to be bounded. This means that the second term of the right-hand side of eq. 1)2.6)1 vanishes 
as e — > +0. Hence we have 



lim ReW(e) = x x = ReUi. 



(2.7) 



Next, we shall consider the imaginary part of W(e), Im W(e). Because of the relation 



arg l e Ul/£ + e U2/e 



arctan 



e xi/e si n (y 1 / e ) + e X2//£ sin(j/2/e) 



e x 1 /e cos(yi/e) + e Z2 / £ cos(?/2/e) 
we have for an infinitesimal value of e and x\ > x 2 , 

,Ui/e + Q U 2 /e 



(2. 



arg l e 



2/1 2/1 
~ arctan tan — = — . 



s £ 

Multiplying e to the both sides of 1)2.9)1 and taking the limit e — > +0, we find 



eargle^^ + e^ 2 ^ 



?/i = Im U\, as e — > +0. 



Hereby from eqs. 1)2.7)1 and 1)2.10)1 . we obtain 

lim W(e) = U\, for Re Ut > Re U 2 . 
By similar procedures, we can see the following holds: 



(2.9) 



(2.10) 



(2.11a) 



lim W{e) = U 2 , for Re Ui <ReU 2 . 



(2.11b) 
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Lastly, let us consider the formula for the case x\ = Xi. For this case, A defined in (|2,5j) 
is given by 

a yi - v% 

A = cos . 

As long as A ^ — 1, that is, a condition 

e + ,f ~ Y\ » ^integer, (2.12) 
(2n + 1J7T 

holds, we can see that (1 + A) £ / 2 goes to unity under e — > 0, and the second term of eq. (|2,6|) 
vanishes in the same limit. Then, assuming ()2.12j) in calculating e — > +0, Re W(e) yields the 
same result as in the case cci 7^ x 2 . As for ImW(e), eq. (|2.8|) reduces to 

arg (e^ + ^\ = arctan [ ^1^+^/^ 
V J Lcos(yi/e) + cos{y 2 /£)_ 

, , VX+V2 VX+V2 , 1 „s 

= arctan tan — = — , (2.16) 



and we arrive at a result 



lim W{e) = Ul + U2 : for ReC/i = ReC/ 2 . (2.14) 
e^+o 2 



Collecting the results (|2,11|) and (|2.14|) . we obtain 

Ui (ReC/x > Ref/ 2 ), 



lim eloeu 



l/ 2 (ReUi < ReU 2 ), (2.15) 

-^4— ^ (ReC/i =Ref/ 2 ). 

This relation is considered to be a generalization of max relation defined for real variables. 
Hereafter, we shall write this operation as A4(f7i,{7 2 ) to distinguish from max(a,6) for real 
numbers a and b. 

3. Operation Rules of M. and Its Generalization to More Than Two Variables 

We shall cmopare the operation rules of M with the max operations for real variables. 
Hereafter in this section, we assume that A, B and C are complex numbers. First, we can 
immediately see from the definition, the commutative law, 

M(A,B)=M(B,A), (3.1) 

is satisfied as is the case for max operation for real numbers. Secondly, by taking logarithm 
of a relation 

a(b + c) = ab + ac, 
and setting e — * +0, we can see that the distributive law 

M(A + B,A + C) = A + M(B,C), (3.2) 
is valid for arbitrary variables. 
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However, the associative law, 

M{A,M(B,C)) = M(M(A,B),C), 
is true at almost all values except at the case 

Re A = ReB = ReC. 
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(3.3a) 



(3.3b) 



This is because the result of taking averages usually depends on the order of operations. This 
situation leads to the difficulty to construct the max-plus algebra for three or more complex 
variables. Hereafter in this section, we are going to define this operation on general number 
of variables. 

Let us consider a summation of ./V variables u±, . . . , u^: 

U = til + • • • + UN- 

As in the previous section, we set the variables as 

Uj = e Uj ' £ , Uj = Xj + it/j, 

(3.4) 

j = l,...,N. 

If the number of the terms which have the largest real parts is two or less, things are the same 
as the previous section, and we have 

M(U U ...,U N ) 

Uk, Uk has the largest real part, (3-5) 
U k + Ui 

, Uk and U\ have the largest real parts. 

We assume that U±, . . . , Uk (K > 3) have the same largest real parts, denoted as x. Calculating 
e log u under an assumption that e is infinitesimal, we have 

elogu = e log (ui H h un) 

~ e In 



e Ui/e + ... + e U K /e 



(3.6) 



is arg I e 



+ 



+ e 



U K /e 



because all the terms e Uj ^ £ (K + 1 < j < N) can be neglected due to the fact 



0, as e — > +0. 



Similarly in the previous section, we find an approximation for infinitesimal e as 

K 



Re (e log u) 



e 



lne 2a;/e + ln 



I cos — + i sm — 

5=1 



x + ^lnK+ |ln(l + A), 



(3.7a) 
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where 

A.lf cos^. (3.7b) 

K ^ £ 

3,1=1 

We can see that the second term of the right-hand side of (|3,7aj) goes to zero under e — > 0, 
and it is also the case for the third term, if e satisfies 

lf>s^/-l. (3.8) 

U=l 

As for the imaginary part of (|3.6|) . because Ui, . . . , Uk have the same real part, we have 

Im (e log it) ~ e arg(e c/l/e H h e Uli/e ) 



£&rg(e iyi/£ + --- + e iyK/£ ) 



Therefore, we obtain a relation 



e log u — > x + i lim e arg ( e iyi ^ £ + • • • + e iyK ^ £ 

From the discussions above, we have derived the definition of A4 for complex variables 
Ui, . . . , Un- We assume that K of them, Uj x , . . . , Uj K , have the same largest real part whose 
value is x, 



U jk =x + iy jk , k = 1, . . . , K. (3.9a) 



Then we have 



M(U 1 ,...,U N ) 



( 



U h , for if = 1, 

/K \ (3.9b) 

x + i lim e arg 2%W e ], forK>2. 

This is the operation of Ai for the general number of complex variables. This definition 
is consistent with eq. ()3.5|) . Unfortunately, although we have defined the hA operation, the 
associative law (j3.3a|) is not satisfied. But if the imaginary parts of Ua (1 < j < K) is 
distributed in equal distance, the latter case of ()3.9b|) can be reduced to 

M(Ux, U K ) = ^p{U\-\ 1- U K ) ■ (3.9c) 

The proof of H3.9cj) is given in the Appendix. 

4. Application of the Ai operation to the Formula of DFT 

In this section, as an application, we shall apply M. operation for DFT. We shall consider 
a variable Ukj, which is a discretization of a function u(x,t), on a region. 

—N < k < N, -M <1<M. (4.1) 
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If we apply DFT on Ukj, we have 

E^-k ( 2irikn 2irilm \ 

^ ^ wexp UivTi + 2MTiJ' (42) 

n=-Nm=-M v 7 

where c„ jm is the Fourier coefficient. 
Now we define 

2nkn 2-Klm 

+ — — -, (4.3a) 



2N + 1 2M + 1 ' 
and e as 

e = J^. (4.3b) 
2iV + 1 v ; 

In order to transform the phase of DFT in the complex max-plus transformation, we introduce 
new coordinates, defined by 

£ = en, v = ek, 

(4.4) 



(4.5) 



(4.6) 



2N + 1 , 2N + 1 

77 = em\ , t = el\ . 

' V 2M + 1 V 2M + 1 

We call these variables stretched coordinates. Substituting (|4.4f) into (|4.3a|) . we have 

= & + T]T 

e 

When we apply the max-plus transformation to u^i and c„ im as 

u k ,i = e u ^ /£ , Cn tm = e c «." /£ , 

U V>T , C^ v : complex numbers, 

and setting e — > +0, we can get a max-plus equation from eq. 1)4. 2JI . The limit e — > +0 is 
equivalent to N — » oo under a condition where M/N is finite. We shall assume that the 
stretched coordinates remain to be constant as N — > oo. As we can see from eqs. 1)4.4)1 . the 
stretched coordinates have discrete values numbered with integers. Then applying suitable 
scale transformations for these variables, we can assume that the values of £ and r\ are restricted 
to be integers. If we express C^ :V = x^ :V + iy^ :V , we have 

Applying the max-plus transformation for complex variables defined in the previous sections, 
we obtain 

U V)T = M ({x^ v + i0£,n )IV r}) , 

(4.8) 

<l>Z,ri,v,T = 2/£,jj +& + f]T. 

Equation ()4.8)1 means that the application of max-plus transformation to DFT yields the 
complex number which have the largest amplitude of discrete Fourier coefficient. We note 
that the trajectory on ^r-plane, U V)T = const., is presented by the line <j>^,rj, v ,r = const., for 
prescribed £ and 77. 
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5. Concluding Remarks and Discussions 

In this paper, we have proposed a method to construct a max-plus equation for complex 
variables. By calculating logarithm of variables explicitly for complex variables, a novel opera- 
tion Ai of max-plus type has been introduced. The conventional max-plus transformation has 
a difficulty that we cannot apply max-plus transformation to non-positive definite vari- 
ables and subtraction terms. For example, an ultradiscretization of the sine-Gordon equation 
has tried by applying max-plus transformation to a discrete analogue of sine-Gordon equation 
to avoid this problem. 6 But it seems that this result owes much to the good structure of 
the original equation. We are convinced that the operation given in this paper dissolves such 
difficulty and enables us to generate discrete systems from difference equations automatically. 

Recently, a way of generalization of max algebra which allows subtraction operation on 
general field was introduced. 7 The relation between the results listed in this reference and the 
one presented in this paper is not clear, but the comparison of these methods is one of the 
worthwhile problems. 

Furthermore, we have applied Ai to DFT. Introducing stretched coordinates clarified the 
relation between DFT and complex max-plus transformation. Since discrete Fourier transfor- 
mation is a common method to solve differential equations, we can expect that solutions of 
various difference equations are easily discretized by virtue of this method. Moreover, we can 
expect that Ai operation enables us to discretize other types of integral transformation, such 
as Laplace transformation, and this will yield a systematic way to solve full discrete systems. 

At the end of this paper, we shall point out several problems concerning the results. The 
definition of the operation Ai for general number of variables (|3.9f) does not allow assosiative 
law, except for the special case shown in eq. (|3.9c[l . This relation still includes limiting process, 
and it seems to be difficult to derive the result of transformation for variables given arbitrary. 
The construction of algebras which satisfy accociative law is a future problem. 
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Appendix: 

We shall prove the formula 



M(U U ... 



U K ) 



Ui + --- + U k 
K 
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under the condition where all the Uj's have the same real part and their imaginary parts are 
distributed in equal distance. Let us set 



Uj = x + iy j} j = l,...,K, (A-2a) 

where 

yj = a + P(j-l). (A-2b) 
The summation of (|3.9b|) is that of a geometric series with common ratio e 1 ^ . Then we have 

i=i 3=1 
sm(/3K/2e) 



K K 



sin(/3/2e 

The argument of this quantity is derived as 

arg (j^wA 



a + 



1 - eW £ 
P(K-l) 



P(K - 1) 



Under the condition (IA-2[) . we find from eq. (|A-4|) that 



M(U U ...,U K ) =x + i 



a + 



P(K-1) 



(A-3) 



(A-4) 



(A-5) 



This is no other than the arithmetic mean of Ux, . . . , Uk defined in (|A-2|) . Hereby we have 
proved eq. (|3.9c|) . 
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